Estrogen receptor 1 chromatin profiling in human breast tumors reveals high inter-patient heterogeneity with enrichment of risk SNPs and enhancer activity at most-conserved regions.

Estrogen Receptor 1 (ESR1; also known as ERα, encoded by ESR1 gene) is the main driver and prime drug target in luminal breast cancer. ESR1 chromatin binding is extensively studied in cell lines and a limited number of human tumors, using consensi of peaks shared among samples. However, little is known about inter-tumor heterogeneity of ESR1 chromatin action, along with its biological implications. Here, we use a large set of ESR1 ChIP-seq data from 70 ESR1+ breast cancers to explore inter-patient heterogeneity in ESR1 DNA binding to reveal a striking inter-tumor heterogeneity of ESR1 action. Of note, commonly shared ESR1 sites show the highest estrogen-driven enhancer activity and are most engaged in long-range chromatin interactions. In addition, the most commonly shared ESR1-occupied enhancers are enriched for breast cancer risk SNP loci. We experimentally confirm SNVs to impact chromatin binding potential for ESR1 and its pioneer factor FOXA1. Finally, in the TCGA breast cancer cohort, we can confirm these variations to associate with differences in expression for the target gene. Cumulatively, we reveal a natural hierarchy of ESR1-chromatin interactions in breast cancers within a highly heterogeneous inter-tumor ESR1 landscape, with the most common shared regions being most active and affected by germline functional risk SNPs for breast cancer development.

GFI1B and LSD1 repress myeloid traits during megakaryocyte differentiation.

The transcription factor Growth Factor Independence 1B (GFI1B) recruits Lysine Specific Demethylase 1 A (LSD1/KDM1A) to stimulate gene programs relevant for megakaryocyte and platelet biology. Inherited pathogenic GFI1B variants result in thrombocytopenia and bleeding propensities with varying intensity. Whether these affect similar gene programs is unknow. Here we studied transcriptomic effects of four patient-derived GFI1B variants (GFI1BT174N,H181Y,R184P,Q287*) in MEG01 megakaryoblasts. Compared to normal GFI1B, each variant affected different gene programs with GFI1BQ287* uniquely failing to repress myeloid traits. In line with this, single cell RNA-sequencing of induced pluripotent stem cell (iPSC)-derived megakaryocytes revealed a 4.5-fold decrease in the megakaryocyte/myeloid cell ratio in GFI1BQ287* versus normal conditions. Inhibiting the GFI1B-LSD1 interaction with small molecule GSK-LSD1 resulted in activation of myeloid genes in normal iPSC-derived megakaryocytes similar to what was observed for GFI1BQ287* iPSC-derived megakaryocytes. Thus, GFI1B and LSD1 facilitate gene programs relevant for megakaryopoiesis while simultaneously repressing programs that induce myeloid differentiation.

Decoding chromatin states by proteomic profiling of nucleosome readers.

DNA and histone modifications combine into characteristic patterns that demarcate functional regions of the genome1,2. While many 'readers' of individual modifications have been described3-5, how chromatin states comprising composite modification signatures, histone variants and internucleosomal linker DNA are interpreted is a major open question. Here we use a multidimensional proteomics strategy to systematically examine the interaction of around 2,000 nuclear proteins with over 80 modified dinucleosomes representing promoter, enhancer and heterochromatin states. By deconvoluting complex nucleosome-binding profiles into networks of co-regulated proteins and distinct nucleosomal features driving protein recruitment or exclusion, we show comprehensively how chromatin states are decoded by chromatin readers. We find highly distinctive binding responses to different features, many factors that recognize multiple features, and that nucleosomal modifications and linker DNA operate largely independently in regulating protein binding to chromatin. Our online resource, the Modification Atlas of Regulation by Chromatin States (MARCS), provides in-depth analysis tools to engage with our results and advance the discovery of fundamental principles of genome regulation by chromatin states.

Inhibition of CSF-1R and IL-6R prevents conversion of cDC2s into immune incompetent tumor-induced DC3s boosting DC-driven therapy potential.

The human dendritic cell (DC) family has recently been expanded by CD1c+CD14+CD163+ DCs, introduced as DC3s. DC3s are found in tumors and peripheral blood of cancer patients. Here, we report elevated frequencies of CD14+ cDC2s, which restore to normal frequencies after tumor resection, in non-small cell lung cancer patients. These CD14+ cDC2s phenotypically resemble DC3s and exhibit increased PD-L1, MERTK, IL-10, and IDO expression, consistent with inferior T cell activation ability compared with CD14- cDC2s. In melanoma patients undergoing CD1c+ DC vaccinations, increased CD1c+CD14+ DC frequencies correlate with reduced survival. We demonstrate conversion of CD5+/-CD1c+CD14- cDC2s to CD14+ cDC2s by tumor-associated factors, whereas monocytes failed to express CD1c under similar conditions. Targeted proteomics identified IL-6 and M-CSF as dominant drivers, and we show that IL-6R and CSF1R inhibition prevents tumor-induced CD14+ cDC2s. Together, this indicates cDC2s as direct pre-cursors of DC3-like CD1c+CD14+ DCs and provides insights into the importance and modulation of CD14+ DC3s in anti-tumor immune responses.

A chromatin-regulated biphasic circuit coordinates IL-1ß-mediated inflammation.

Inflammation is characterized by a biphasic cycle consisting initially of a proinflammatory phase that is subsequently resolved by anti-inflammatory processes. Interleukin-1ß (IL-1ß) is a master regulator of proinflammation and is encoded within the same topologically associating domain (TAD) as IL-37, which is an anti-inflammatory cytokine that opposes the function of IL-1ß. Within this TAD, we identified a long noncoding RNA called AMANZI, which negatively regulates IL-1ß expression and trained immunity through the induction of IL37 transcription. We found that the activation of IL37 occurs through the formation of a dynamic long-range chromatin contact that leads to the temporal delay of anti-inflammatory responses. The common variant rs16944 present in AMANZI augments this regulatory circuit, predisposing individuals to enhanced proinflammation or immunosuppression. Our work illuminates a chromatin-mediated biphasic circuit coordinating expression of IL-1ß and IL-37, thereby regulating two functionally opposed states of inflammation from within a single TAD.

Acid ceramidase regulates innate immune memory.

Innate immune memory, also called "trained immunity," is a functional state of myeloid cells enabling enhanced immune responses. This phenomenon is important for host defense, but also plays a role in various immune-mediated conditions. We show that exogenously administered sphingolipids and inhibition of sphingolipid metabolizing enzymes modulate trained immunity. In particular, we reveal that acid ceramidase, an enzyme that converts ceramide to sphingosine, is a potent regulator of trained immunity. We show that acid ceramidase regulates the transcription of histone-modifying enzymes, resulting in profound changes in histone 3 lysine 27 acetylation and histone 3 lysine 4 trimethylation. We confirm our findings by identifying single-nucleotide polymorphisms in the region of ASAH1, the gene encoding acid ceramidase, that are associated with the trained immunity cytokine response. Our findings reveal an immunomodulatory effect of sphingolipids and identify acid ceramidase as a relevant therapeutic target to modulate trained immunity responses in innate immune-driven disorders.

Breast cancer risk SNPs converge on estrogen receptor binding sites commonly shared between breast tumors to locally alter estrogen signalling output.

Estrogen Receptor alpha (ERα) is the main driver and prime drug target in luminal breast. ERα chromatin binding is extensively studied in cell lines and a limited number of human tumors, using consensi of peaks shared among samples. However, little is known about inter-tumor heterogeneity of ERα chromatin action, along with its biological implications. Here, we use a large set of ERα ChIP-seq data from 70 ERα+ breast cancers to explore inter-patient heterogeneity in ERα DNA binding, to reveal a striking inter-tumor heterogeneity of ERα action. Interestingly, commonly-shared ERα sites showed the highest estrogen-driven enhancer activity and were most-engaged in long-range chromatin interactions. In addition, the most-commonly shared ERα-occupied enhancers were enriched for breast cancer risk SNP loci. We experimentally confirm SNVs to impact chromatin binding potential for ERα and its pioneer factor FOXA1. Finally, in the TCGA breast cancer cohort, we could confirm these variations to associate with differences in expression for the target gene. Cumulatively, we reveal a natural hierarchy of ERα-chromatin interactions in breast cancers within a highly heterogeneous inter-tumor ERα landscape, with the most-common shared regions being most active and affected by germline functional risk SNPs for breast cancer development.

DDA1, a novel factor in transcription-coupled repair, modulates CRL4CSA dynamics at DNA damage-stalled RNA polymerase II.

Transcription-blocking DNA lesions are specifically targeted by transcription-coupled nucleotide excision repair (TC-NER), which removes a broad spectrum of DNA lesions to preserve transcriptional output and thereby cellular homeostasis to counteract aging. TC-NER is initiated by the stalling of RNA polymerase II at DNA lesions, which triggers the assembly of the TC-NER-specific proteins CSA, CSB and UVSSA. CSA, a WD40-repeat containing protein, is the substrate receptor subunit of a cullin-RING ubiquitin ligase complex composed of DDB1, CUL4A/B and RBX1 (CRL4CSA). Although ubiquitination of several TC-NER proteins by CRL4CSA has been reported, it is still unknown how this complex is regulated. To unravel the dynamic molecular interactions and the regulation of this complex, we applied a single-step protein-complex isolation coupled to mass spectrometry analysis and identified DDA1 as a CSA interacting protein. Cryo-EM analysis showed that DDA1 is an integral component of the CRL4CSA complex. Functional analysis revealed that DDA1 coordinates ubiquitination dynamics during TC-NER and is required for efficient turnover and progression of this process.

A cyclin-dependent kinase-mediated phosphorylation switch of disordered protein condensation.

Cell cycle transitions result from global changes in protein phosphorylation states triggered by cyclin-dependent kinases (CDKs). To understand how this complexity produces an ordered and rapid cellular reorganisation, we generated a high-resolution map of changing phosphosites throughout unperturbed early cell cycles in single Xenopus embryos, derived the emergent principles through systems biology analysis, and tested them by biophysical modelling and biochemical experiments. We found that most dynamic phosphosites share two key characteristics: they occur on highly disordered proteins that localise to membraneless organelles, and are CDK targets. Furthermore, CDK-mediated multisite phosphorylation can switch homotypic interactions of such proteins between favourable and inhibitory modes for biomolecular condensate formation. These results provide insight into the molecular mechanisms and kinetics of mitotic cellular reorganisation.

Transcriptomics and proteomics reveal distinct biology for lymph node metastases and tumour deposits in colorectal cancer.

Both lymph node metastases (LNMs) and tumour deposits (TDs) are included in colorectal cancer (CRC) staging, although knowledge regarding their biological background is lacking. This study aimed to compare the biology of these prognostic features, which is essential for a better understanding of their role in CRC spread. Spatially resolved transcriptomic analysis using digital spatial profiling was performed on TDs and LNMs from 10 CRC patients using 1,388 RNA targets, for the tumour cells and tumour microenvironment. Shotgun proteomics identified 5,578 proteins in 12 different patients. Differences in RNA and protein expression were analysed, and spatial deconvolution was performed. Image-based consensus molecular subtype (imCMS) analysis was performed on all TDs and LNMs included in the study. Transcriptome and proteome profiles identified distinct clusters for TDs and LNMs in both the tumour and tumour microenvironment segment, with upregulation of matrix remodelling, cell adhesion/motility, and epithelial-mesenchymal transition (EMT) in TDs (all p < 0.05). Spatial deconvolution showed a significantly increased abundance of fibroblasts, macrophages, and regulatory T-cells (p < 0.05) in TDs. Consistent with a higher fibroblast and EMT component, imCMS classified 62% of TDs as poor prognosis subtype CMS4 compared to 36% of LNMs (p < 0.05). Compared to LNMs, TDs have a more invasive state involving a distinct tumour microenvironment and upregulation of EMT, which are reflected in a more frequent histological classification of TDs as CMS4. These results emphasise the heterogeneity of locoregional spread and the fact that TDs should merit more attention both in future research and during staging. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Galectin-9 interacts with Vamp-3 to regulate cytokine secretion in dendritic cells.

Intracellular vesicle transport is essential for cellular homeostasis and is partially mediated by SNARE proteins. Endosomal trafficking to the plasma membrane ensures cytokine secretion in dendritic cells (DCs) and the initiation of immune responses. Despite its critical importance, the specific molecular components that regulate DC cytokine secretion are poorly characterised. Galectin-9, a ß-galactoside-binding protein, has emerged as a novel cellular modulator although its exact intracellular roles in regulating (immune) cell homeostasis and vesicle transport are virtually unknown. We investigated galectin-9 function in primary human DCs and report that galectin-9 is essential for intracellular cytokine trafficking to the cell surface. Galectin-9-depleted DCs accumulate cytokine-containing vesicles in the Golgi complex that eventually undergo lysosomal degradation. We observed galectin-9 to molecularly interact with Vamp-3 using immunoprecipitation-mass-spectrometry and identified galectin-9 was required for rerouting Vamp-3-containing endosomes upon DC activation as the underlying mechanism. Overall, this study identifies galectin-9 as a necessary mechanistic component for intracellular trafficking. This may impact our general understanding of vesicle transport and sheds new light into the multiple roles galectins play in governing cell function.

STPP-UP: An alternative method for drug target identification using protein thermal stability.

Thermal proteome profiling (TPP) has significantly advanced the field of drug discovery by facilitating proteome-wide identification of drug targets and off-targets. However, TPP has not been widely applied for high-throughput drug screenings, since the method is labor intensive and requires a lot of measurement time on a mass spectrometer. Here, we present Single-tube TPP with Uniform Progression (STPP-UP), which significantly reduces both the amount of required input material and measurement time, while retaining the ability to identify drug targets for compounds of interest. By using incremental heating of a single sample, changes in protein thermal stability across a range of temperatures can be assessed, while alleviating the need to measure multiple samples heated to different temperatures. We demonstrate that STPP-UP is able to identify the direct interactors for anticancer drugs in both human and mice cells. In summary, the STPP-UP methodology represents a useful tool to advance drug discovery and drug repurposing efforts.

Live cell transcription-coupled nucleotide excision repair dynamics revisited.

Transcription-blocking lesions are specifically targeted by transcription-coupled nucleotide excision repair (TC-NER), which prevents DNA damage-induced cellular toxicity and maintains proper transcriptional processes. TC-NER is initiated by the stalling of RNA polymerase II (RNAPII), which triggers the assembly of TC-NER-specific proteins, namely CSB, CSA and UVSSA, which collectively control and drive TC-NER progression. Previous research has revealed molecular functions for these proteins, however, exact mechanisms governing the initiation and regulation of TC-NER, particularly at low UV doses have remained elusive, partly due to technical constraints. In this study, we employ knock-in cell lines designed to target the endogenous CSB gene locus with mClover, a GFP variant. Through live cell imaging, we uncover the intricate molecular dynamics of CSB in response to physiologically relevant UV doses. We showed that the DNA damage-induced association of CSB with chromatin is tightly regulated by the CSA-containing ubiquitin-ligase CRL complex (CRL4CSA). Combining the CSB-mClover knock-in cell line with SILAC-based GFP-mediated complex isolation and mass-spectrometry-based proteomics, revealed novel putative CSB interactors as well as discernible variations in complex composition during distinct stages of TC-NER progression. Our work not only provides molecular insight into TC-NER, but also illustrates the versatility of endogenously tagging fluorescent and affinity tags.

Myeloid-derived suppressor cells and tolerogenic dendritic cells are distinctively induced by PI3K and Wnt signaling pathways.

Imbalanced immune responses are a prominent hallmark of cancer and autoimmunity. Myeloid cells can be overly suppressive, inhibiting protective immune responses or inactive not controlling autoreactive immune cells. Understanding the mechanisms that induce suppressive myeloid cells, such as myeloid-derived suppressor cells (MDSCs) and tolerogenic dendritic cells (TolDCs), can facilitate the development of immune-restoring therapeutic approaches. MDSCs are a major barrier for effective cancer immunotherapy by suppressing antitumor immune responses in cancer patients. TolDCs are administered to patients to promote immune tolerance with the intent to control autoimmune disease. Here, we investigated the development and suppressive/tolerogenic activity of human MDSCs and TolDCs to gain insight into signaling pathways that drive immunosuppression in these different myeloid subsets. Moreover, monocyte-derived MDSCs (M-MDSCs) generated in vitro were compared to M-MDSCs isolated from head-and-neck squamous cell carcinoma patients. PI3K-AKT signaling was identified as being crucial for the induction of human M-MDSCs. PI3K inhibition prevented the downregulation of HLA-DR and the upregulation of reactive oxygen species and MerTK. In addition, we show that the suppressive activity of dexamethasone-induced TolDCs is induced by ß-catenin-dependent Wnt signaling. The identification of PI3K-AKT and Wnt signal transduction pathways as respective inducers of the immunomodulatory capacity of M-MDSCs and TolDCs provides opportunities to overcome suppressive myeloid cells in cancer patients and optimize therapeutic application of TolDCs. Lastly, the observed similarities between generated- and patient-derived M-MDSCs support the use of in vitro-generated M-MDSCs as powerful model to investigate the functionality of human MDSCs.

A novel antifolate suppresses growth of FPGS-deficient cells and overcomes methotrexate resistance.

Cancer cells make extensive use of the folate cycle to sustain increased anabolic metabolism. Multiple chemotherapeutic drugs interfere with the folate cycle, including methotrexate and 5-fluorouracil that are commonly applied for the treatment of leukemia and colorectal cancer (CRC), respectively. Despite high success rates, therapy-induced resistance causes relapse at later disease stages. Depletion of folylpolyglutamate synthetase (FPGS), which normally promotes intracellular accumulation and activity of natural folates and methotrexate, is linked to methotrexate and 5-fluorouracil resistance and its association with relapse illustrates the need for improved intervention strategies. Here, we describe a novel antifolate (C1) that, like methotrexate, potently inhibits dihydrofolate reductase and downstream one-carbon metabolism. Contrary to methotrexate, C1 displays optimal efficacy in FPGS-deficient contexts, due to decreased competition with intracellular folates for interaction with dihydrofolate reductase. We show that FPGS-deficient patient-derived CRC organoids display enhanced sensitivity to C1, whereas FPGS-high CRC organoids are more sensitive to methotrexate. Our results argue that polyglutamylation-independent antifolates can be applied to exert selective pressure on FPGS-deficient cells during chemotherapy, using a vulnerability created by polyglutamylation deficiency.

Determining DNA-Protein Binding Affinities and Specificities from Crude Lysates Using a Combined SILAC/TMT Labeling Strategy.

In recent years, quantitative mass spectrometry-based interaction proteomics technology has proven very useful in identifying specific DNA-protein interactions using single pull-downs from crude lysates. Here, we applied a SILAC/TMT-based higher-order multiplexing approach to develop an interaction proteomics workflow called Protein-nucleic acid Affinity and Specificity quantification by MAss spectrometry in Nuclear extracts or PASMAN. In PASMAN, DNA pull-downs using a concentration range of specific and control DNA baits are performed in SILAC-labeled nuclear extracts. MS1-based quantification to determine specific DNA-protein interactions is then combined with sequential TMT-based quantification of fragmented SILAC peptides, allowing the generation of Hill-like curves and determination of apparent binding affinities. We benchmarked PASMAN using the SP/KLF motif and further applied it to gain insights into two CGCG-containing consensus DNA motifs. These motifs are recognized by two BEN domain-containing proteins, BANP and BEND3, which we find to interact with these motifs with distinct affinities. Finally, we profiled the BEND3 proximal proteome, revealing the NuRD complex as the major BEND3 proximal protein complex in vivo. In summary, PASMAN represents, to our knowledge, the first higher-order multiplexing-based interaction proteomics method that can be used to decipher specific DNA-protein interactions and their apparent affinities in various biological and pathological contexts.

THRONCAT: metabolic labeling of newly synthesized proteins using a bioorthogonal threonine analog.

Profiling the nascent cellular proteome and capturing early proteomic changes in response to external stimuli provides valuable insights into cellular physiology. Existing metabolic protein labeling approaches based on bioorthogonal methionine- or puromycin analogs allow for the selective visualization and enrichment of newly synthesized proteins. However, their applications are limited as they often require methionine-free conditions, auxotrophic cells and/or are toxic to cells. Here, we introduce THRONCAT, a threonine-derived non-canonical amino acid tagging method based on the bioorthogonal threonine analog ß-ethynylserine (ßES) that enables efficient labeling of the nascent proteome in complete growth media within minutes. We use THRONCAT for the visualization and enrichment of nascent proteins in bacteria, mammalian cells and Drosophila melanogaster. We profile immediate proteome dynamics of B-cells in response to B-cell receptor activation simply by adding ßES to the culture medium, demonstrating the ease-of-use of the method and its potential to address diverse biological questions. In addition, using a Drosophila model of Charcot-Marie-Tooth peripheral neuropathy, we show that THRONCAT enables visualization and quantification of relative protein synthesis rates in specific cell types in vivo.

Dynamics and Composition of Small Heat Shock Protein Condensates and Aggregates.

Small heat shock proteins (sHSPs) are essential ATP-independent chaperones that protect the cellular proteome. These proteins assemble into polydisperse oligomeric structures, the composition of which dramatically affects their chaperone activity. The biomolecular consequences of variations in sHSP ratios, especially inside living cells, remain elusive. Here, we study the consequences of altering the relative expression levels of HspB2 and HspB3 in HEK293T cells. These chaperones are partners in a hetero-oligomeric complex, and genetic mutations that abolish their mutual interaction are associated with myopathic disorders. HspB2 displays three distinct phenotypes when co-expressed with HspB3 at varying ratios. Expression of HspB2 alone leads to formation of liquid nuclear condensates, while shifting the stoichiometry towards HspB3 resulted in the formation of large solid-like aggregates. Only cells co-expressing HspB2 with a limited amount of HspB3 formed fully soluble complexes that were distributed homogeneously throughout the nucleus. Strikingly, both condensates and aggregates were reversible, as shifting the HspB2:HspB3 balance in situ resulted in dissolution of these structures. To uncover the molecular composition of HspB2 condensates and aggregates, we used APEX-mediated proximity labelling. Most proteins interact transiently with the condensates and were neither enriched nor depleted in these cells. In contrast, we found that HspB2:HspB3 aggregates sequestered several disordered proteins and autophagy factors, suggesting that the cell is actively attempting to clear these aggregates. This study presents a striking example of how changes in the relative expression levels of interacting proteins affects their phase behavior. Our approach could be applied to study the role of protein stoichiometry and the influence of client binding on phase behavior in other biomolecular condensates and aggregates.

TGF-ß Type I Receptor Signaling in Melanoma Liver Metastases Increases Metastatic Outgrowth.

Despite advances in treatment for metastatic melanoma patients, patients with liver metastasis have an unfavorable prognosis. A better understanding of the development of liver metastasis is needed. The multifunctional cytokine Transforming Growth Factor ß (TGF-ß) plays various roles in melanoma tumors and metastasis, affecting both tumor cells and cells from the surrounding tumor microenvironment. To study the role of TGF-ß in melanoma liver metastasis, we created a model to activate or repress the TGF-ß receptor pathway in vitro and in vivo in an inducible manner. For this, we engineered B16F10 melanoma cells to have inducible ectopic expression of a constitutively active (ca) or kinase-inactive (ki) TGF-ß receptor I, also termed activin receptor-like kinase (ALK5). In vitro, stimulation with TGF-ß signaling and ectopic caALK5 expression reduced B16F10 cell proliferation and migration. Contrasting results were found in vivo; sustained caALK5 expression in B16F10 cells in vivo increased the metastatic outgrowth in liver. Blocking microenvironmental TGF-ß did not affect metastatic liver outgrowth of both control and caALK5 expressing B16F10 cells. Upon characterizing the tumor microenvironment of control and caALk5 expressing B16F10 tumors, we observed reduced (cytotoxic) T cell presence and infiltration, as well as an increase in bone marrow-derived macrophages in caALK5 expressing B16F10 tumors. This suggests that caALK5 expression in B16F10 cells induces changes in the tumor microenvironment. A comparison of newly synthesized secreted proteins upon caALK5 expression by B16F10 cells revealed increased secretion of matrix remodeling proteins. Our results show that TGF-ß receptor activation in B16F10 melanoma cells can increase metastatic outgrowth in liver in vivo, possibly through remodeling of the tumor microenvironment leading to altered infiltration of immune cells. These results provide insights in the role of TGF-ß signaling in B16F10 liver metastasis and could have implications regarding the use of TGF-ß inhibitors for the treatment of melanoma patients with liver metastasis.